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Abstract 

Iron nitride Fe4N is studied by full-potential LAPW method. Structure 
parameters, electronic and magnetic properties as well as hyperfine interaction 
parameters are obtained. We observe perfect agreement with experimental re- 
sults. Hypothetical Fe4N structure was also calculated to study the influence of 
disordering effects on parameters of Mossbauer spectra. We performed detailed 
analysis of EFG formation on Fe nuclei including magnetization effects. We show 
that the formation of N-Fe-N local configuration is energetically favourable in 
nitrogen austenites. 



1 Introduction 

Face centered cubic (fee) iron-based alloys are widely used for developing of stainless 
austenitic steels especially for using in critical temperature ranges, aggressive environ- 
ment and other severe external conditions. Doping of these steels with light interstitial 
impurities (C,N) influence mechanics and kinetics of structure phase transitions in Fe- 
based alloys. Distribution of carbon and nitrogen atoms in solid solutions influence 
electrical and mechanical properties of alloys. Nitrogen doping enables to solve the 
problem of the strengthening of stainless steels. 

Investigation of the influence of nitrogen on physical properties of multicomponent 
systems is a complicated problem. The solution of this problem should be made in 
several stages. On the first stage it seems important to study electronic structure of 
iron nitride Fe4N, because binary 7-FeNa. alloy can be modeled by non-stoichiometric 
nitride 7-Fe4Ni_a.. 

There are a lot of experimental data about atomic and electronic structure of Fe4N. 
We believe that Mossbauer spectroscopy gives most interesting information about im- 
purity distribution, electronic structure and magnetic interactions in alloys. Studying 
of hyperfine structure of the energy spectra of nuclei is a powerful tool for investiga- 
tion of interactions of atomic nucleus with local electric and magnetic fields. These 
interactions cause shifts and splittings of nuclear energy levels and enable us to ob- 
tain information about symmetry of charge distribution around the nucleus, about 
electronic configurations of atoms and ions, as well as about peculiarities of atomic 
structure of solids. 

A number of experimental papers show substantial differences in Mossbauer spec- 
tra of binary Fe-N and Fe-C alloys. These differences are believed to be connected 



with different C and N distribution in Fe fee alloys |]T|. In this paper we present cal- 
culation results of hyperfine interaction parameters for iron nitride Fe4N as well as for 
hypothetical Fe4N structure with another distribution of nitrogen atoms. This allows 
us to determine changes in Mossbauer spectra caused by redistribution on nitrogen 
atoms. 

2 Calculation details 

WIEN97 programme package p[, employing Full-potential Linearized Augmented Plane 
Wave (FLAPW) method was used for calculations. As far as FLAPW is an all-electron 
method (unlike pseudopotential methods), it allows to perform calculations of hyper- 
fine interaction parameters from first principles. Obtained theoretical parameters of 
interaction of a nucleus with electric and magnetic fields can be successfully compared 
with parameters of experimental Mossbauer spectra. 

Generalized gradient approximation (GGA) according to Perdew-Burke-Ernzerhof 
model was used for exchange- correlation potential. The radii of atomic spheres 
were chosen as 1.9 a.u. and 1.68 a.u for Fe and N atoms respectively. The accuracy of 
calculation results depends on several basic parameters: number of fC-points in Bril- 
louin zone, number of LM-components and Fourier coefficients in charge density and 
potential decomposition and number of plane waves in interstitial region. The choice 
of the values of these parameters was based on convergence condition. Convergence 
tests gave the value Rmin x Kmax = 8.8, which corresponds to 205 plane waves per 
atom in the basis set. Inside atomic spheres the wave function was decomposed up to 
imax = 12. Charge density and potential was decomposed inside atomic spheres using 
lattice harmonics basis up to L^ax = 6. In the interstitial region Fourier expansion 
was used with 850 coefficients. Calculations were performed for 3000 K-points in the 
Brillouin zone (84 K-points in the irreducible part). The values of all parameters en- 
sure accuracy of 0.1 mRy in total energy of the system. Due to ferromagnetic nature 
of iron nitride all calculations were performed using spin-polarized approximation. 

3 Atomic structure 

The unit cell of iron nitride Fe4N is a unit cell of 7-Fe with a nitrogen atom introduced 
in the centre of the cube (figure [^a). The structure has two symmetry types of Fe 
atoms: Fe2 type forms octahedron around impurity atom and has two impurity atoms 
in the first coordination sphere located at 180°; Fcq type is located in the corners of 
the cell and has no impurity atoms in the first coordination sphere. 

In order to find the value of lattice parameter, corresponding to the minimum total 
energy of the system we performed calculations for five different values of lattice pa- 
rameter. The results were approximated by second-order polynomial using least square 
fit method. Then the value of lattice parameter was obtained analytically. Calculated 
value of lattice parameter C=7.164 a.u. is in good agreement with experimental value 
C=7.17 a.u. [§. Calculations of electronic structure and hyperfine interaction param- 
eters were performed using obtained optimized value of lattice parameter. 

Optimization procedure was also performed for another two systems: 7-Fe (Fe-fcc) 
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Figure 1: The unit cell of iron nitride (a) and of hypothetical Fe4N structure (b) 



and hypothetical Fe4N structure (figure |I|b). Both calculations were also performed in 
spin-polarized approximation. Calculation of 7-Fe seems interesting because it allows 
us to determine the influence of introducing of nitrogen atom on electronic structure 
and hyperfine interaction parameters of Fe-fcc. Optimization of hypothetical Fe4N 
structure (figure |pD) was performed including variation of interatomic Fe-N distance. 
This structure has two types of iron atoms: Fei and Fe2, having one and two N atoms in 
the nearest neighbourhood respectively. We found out that with this type of ordering 
the total energy of the structure is 18 mRy larger, than the total energy of iron nitride. 
This shows that the formation of Fe-N-Fe-N chains (Fe2 configuration) without Fe-N 
pairs (Fei configuration) is energetically favourable. 

From the point of view of atomic structure of Fe-N alloys iron nitride Fe4N is the 
simplest structure where Fe-N chains are combined with Fe-Fe chains. This makes 
possible to perform detailed analysis of the influence of nitrogen on electronic struc- 
ture and hyperfine interaction parameters and to use obtained information later for 
analyzing more complicated Fe-N alloys. This step-by-step approach, based on ap- 
plication of up-to-date ab initio calculation method, appears to be most accurate for 
solving the problem of the influence of impurity atoms on the properties of Fe-based 
alloys. 



Table 1: Cell parameters (a,b,c) and magnetic moments (M) for 7-Fe, iron nitride and model Fe4N 
structure. Cell parameter is presented in atomic units, magnetic moment - in Bohr magnetons. 



Structure 


a, 6, c (a.u) 


a, 5, c (a.u.) 




M{^iB) 




theory 


expt. 


theory 


expt. 


Fe(fcc) 


a = 6.58 


a = 6.8T 


1.04 




Fe4N 


a = 7.164 


a = 7.171" 






Feo 






2.90 


S.O'^ 


Fe2 






2.31 


2.0^= 


Aver. /atom 






2.50 


2.21^= 


Fe4N (model) 


a=7.114 










c = 14.266 








Fei 






2.53 




Fea 






1.98 




Aver. /atom 






2.27 





0", & 
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4 Hyperfine interactions 



Ferromagnetic nature of iron nitride Fe4N makes Mossbauer spectrum of this com- 
pound much more comphcated. The presence of magnetic field leads to combined 
magnetic and electric hyperfine interactions. The Hamiltonian of the system, includ- 
ing these interactions is given by 

?i = -(/ixi?)+e^g,,V^,,, (1) 

where /2 - magnetic moment of the nucleus, H - magnetic field, Qij - tensor of the 
nuclear quadrupole moment, Vij - tensor of electric field gradient (EFG). The princi- 
pal EFG axis and magnetization axis may be non-coUinear, which makes finding of 
Hamiltonian eigenvalues much more complicated. Using eQ/fiH <^ 1 approximation 
and taking into account that 1 = 3/2 for excited state of ^^Fe nucleus, the energy 
levels of the system are given by 




m/ = -/, -/ + 1,. . .,/ - 1,/ 



In this expression m/ - magnetic quantum number, Vzz - the principal component of 
EFG tensor, rj - asymmetry parameter, 6 - angle between the principal axis of EFG 
tensor and the direction of magnetic field H. 

The presence of magnetic field reduces the symmetry of the system. In iron nitride 
Fe atoms, located in the centre of the cube faces are now divided into two types: 
Fe2, for which principal EFG axis is parallel to magnetic field, and Fe2, for which 
this axis is perpendicular to magnetic field (figure |l|a). Mossbauer spectra of these 
two types of Fe atoms will have different quadrupole splitting. So, due to magnetic 
field, we have three different types of iron atoms in iron nitride structure. The total 
Mossbauer spectrum of this compound will consist of eighteen peaks. Twelve of these 
peaks, which belong to Fe2 and Fe'g atoms, will have different shifts due to different 
quadrupole interactions. The ratio of the number of these atomic types in iron nitride 
unit cell is Feo:Fe2:Fe2=l:l:2, which leads to the same ratio of intensities of respective 
Mossbauer peaks. All described features of the spectrum are totally approved by 
experimental data . 




Quadrupole splitting, expressed in mm/s, is given by 

A4mm/s] = 0.5205 x Q[b]Vzz[lQ''\ / m^f-^^^^—^ ^^ ^) ' 

In our calculations we used the value of nuclear quadrupole moment Q^^Fe =0.166. 
This value was determined by Dufek, Blaha and Schwarz [^] by comparing experimental 
quadrupole splitting and theoretical EFG values, calculated by FLAPW method, for 
fourteen different Fe compounds. 

For temperatures above the Curie point (Tc) the system becomes paramagnetic, 
and Fcq atoms will not contribute to Mossbauer spectrum. Fe2 and Fe'g atomic types 
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Table 2: Contributions to Vzzi^O'^^v/m'^) from inside the spheres and from interstitial region. 
Theoretical and experimental values of quadrupole splitting (QS)(mm/s) for Fe atoms with different 
angle 6 between H and principal EFG axis. 



System 


Atom 


Valence 
EFG 


Lattice 
EFG 


Total 
EFG 


V 


9° 


QS 

theory 


QS 

expt. H 


Nitride 




















Fe2 


-3.369 


+0.382 


-2.987 





0° 


-0.25 


-0.26 + 0.02 




Fe^ 


-3.369 


+0.382 


-2.987 





90° 


0.12 


0.12 + 0.02 


Model 




















Fei 


-0.34 


+0.14 


-0.20 





90° 


0.01 






Fea 


-2.41 


+0.62 


-1.79 


0.76 


0° 


-0.16 






Fe^ 


-2.41 


+0.62 


-1.79 


0.76 


90° 


0.08 





become totally identical, and Mossbauer spectrum in this case is a doublet with 
quadrupole splitting 




A^[mm/s] = lMlxQ[b]V,,[10^h/m^]{l + ^] . (4) 



In this paper electric field gradient was calculated according to method, developed 
by Blaha et al [0. According to this approach EFG is calculated on ab-initio basis 



directly from electronic density distribution. We obtain EFG contributions from charge 
density inside atomic spheres {valence EFG) and outside the spheres {lattice EFG). 
Valence EFG is calculated as an integral of the value P2M{f^)/f^ over the atomic sphere 
||11||. The values plm originate from two radial wave functions with I and 



PLM{r)= E j:j:RUr)Ri'm'{r)G^ir, (5) 

E<Ep l,m l',m' 

where Rim{r) - LAPW radial wave functions, G^^™"^' - Gaunt integrals. 

Table H presents calculated principal values of EFG tensor and obtained according 
to (I) values of quadrupole splitting (QS) for Fe2 and Fe'a atoms. The table shows 
that our results are in perfect agreement with experimental data, obtained by Foct 



et al 0] for the sample at T=4K. Excellent agreement with experiment also shows 
that chosen direction of magnetic axis (001) is a direction of real magnetic field in iron 
nitride Fe4N. 

Using (^) we also calculated quadrupole splitting at T > Tc. Our value Ai,=0.50 
mm/s is also in excellent agreement with experimental data of Foct (0.50 mm/s). 



obtained at T=763K 12 



Table ^ also presents calculation results for model Fe4N structure. Magnetization 
axis is assumed to be directed along one of Fe-N-Fe-N chains. Redistribution of ni- 
trogen atoms leads to substantial changes of EFG values. The quadrupole splitting 
of Fe2 atoms is reduced by 1/3 of the previous value when Fei configuration appears 
in the structure. The same effect we observe at T>Tc: quadrupole splitting becomes 
now 0.33 mm/s. We see that quadrupole splitting is very sensitive to redistribution 
of impurity atoms in Fe-fcc. This fact can be used for interpretation of Mossbauer 
spectra of real fcc-alloys with different distribution of interstitial impurity. 
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5 Charge density and electric field gradient 



We performed detailed analysis of EFG formation in iron nitride Fe4N. Table H shows 
that valence component makes 90% contribution to the total value. So, we can focus 
on this component for understanding of EFG origin. 

According to our approach the valent band of the compound (Fe-3s3p3d, N-2s2p) 
was divided to four energy regions according to table |^. These energy intervals are 
formed mostly by different electronic shells of iron and nitrogen atoms. We also ana- 
lyzed EFG contributions from wavefunctions of different symmetry (according to (j^)) 
for each of energy regions. 

Table 3: Principal EFG component (10^^ v/rn^) for iron nitride Fe4N. Contributions to Vzz from 
different energy intervals and from radial wavefunctions of different symmetry. 



Energy 
region 


s — d 


P-P 


d-d 


Total 


Fe-3s3p 


-0.046 


+4.437 


-0.001 


+4.390 


N-2s 


-0.021 


-7.889 


-0.364 


-8.274 


N-2p 


-0.037 


-14.652 


-5.680 


-20.369 


Fe-3d 


+0.092 


+7.892 


+13.121 


+21.105 


Valent 


-0.012 


-10.212 


+7.076 


-3.148 



Table 4: V^z contributions from N-2p and Fe-3ci energy regions for electrons with different spin 
orientation. 



Energy 


Spin 


s — d 


P-P 


d-d 


Total 


region 












N-2p 


up 


-0.015 


-6.535 


-3.021 


-9.571 




down 


-0.022 


-8.117 


-2.659 


-10.798 


Fe-3d 


up 


+0.042 


+2.216 


+ 10.190 


+ 12.448 




down 


+0.050 


+5.676 


+2.931 


+8.657 



First of all let us analyze the electric field gradient, caused by the whole valent 
band (valence EFG). Table |^ shows that s — d contribution is negligibly small for all 
energy regions. EFG is practically formed hj p — p and d — d contributions, which 
have different signs. We see that negative p — p contribution is dominating, although 
it is greatly canceled by positive d — d contribution. 

The analysis of table ^ data allows us to conclude that main contribution to EFG 
formation is made by 2p and 2s nitrogen states, and contribution from 2p states is 
dominant. It should be noted that we are now talking about electrons, which have 
energy of nitrogen 2p states, but spatially are located in Fe2 sphere. Substantial 
asymmetry of spatial distribution of these electrons (figure |]a) gives main negative 
contribution to total EFG. The interesting fact is that more than 70% of charge, 
formed by these electrons, have p— symmetry, and 30% - (i— symmetry relative to the 
centre of Fe2 sphere. This follows from analysis oi p — p and d — d contributions for 
N-2p interval (table The same situation is observed for 2s— electrons of nitrogen, 
but their total EFG contribution is 2.5 times smaller than 2p— electrons contribution. 
The EFG contribution of N-2s— electrons is practically caused only by electrons of 
p— symmetry relative to the centre of Fe2 sphere. 
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Substantial negative EFG, caused by described asymmetry of nitrogen electrons, 
is greatly canceled by positive gradient, caused by own Fe electrons (figure ||b). These 
are first of all Fe-3(i— electrons, which density is distributed in such a way, that greatly 
cancels asymmetry of 2s and 2p electrons of nitrogen. Comparative analysis of p—p and 
d — d contributions shows that more than 60% of asymmetrically distributed electrons 
of Fe-3d energy region have (i— symmetry, and less than 40% - symmetry (table ^ 
relative to the centre of Fe2 sphere. Redistribution of rather localized Fe-3s3p states 
also gives positive EFG contribution, but this contribution is almost 5 times smaller 
than contribution of Sci— states. 

It should be noted that quadrupole splitting was calculated using total EFG value, 
given in table ^, and for calculation of valence EFG contribution summation in (^ was 
done up to Imax = 12. On the other hand, for analysis of partial EFG contributions 
we limited the summation only up to Imax = 2 (due to the absence of "chemical" 
/-states). This is the reason for small differences of valence EFG values in tables |^ 
and 1^. 

The values of local magnetic moments play an important role in EFG formation in 
iron nitride Fe4N. The calculations show that the upper edge of Fe-3d spin down states 
is located 2 eV above the Fermi level. Substantial difference in population of states 
with different spin orientation greatly influence EFG formation. Table ^ gives EFG 
contributions from states with different spin orientation for N-2p and Fe-3d energy 
regions. We notice considerable difference in EFG values, caused by spin up and down 
electrons of Fe-3(i interval. This effect mostly influence the value of d — d contribution 
of this energy region. This contribution is mainly formed by density distribution of 
own Sci— electrons of iron. Spin up contribution is more than 3 times larger than spin 
down contribution for d — d component of Fe-3(i— region. 

The fact that spin down states are partially "pulled" above the Fermi level consid- 
erably decreases compensating positive component, and leads to increasing of absolute 
value of total (negative) EFG. In the absence of this effect we would have considerably 
larger compensating component, which would lead to decreasing of total EFG value 
and, as a consequence, to decreasing of quadrupole splitting in Mossbauer spectra. 




Figure 2: Charge distribution in N-2p (a) and N-2pFe3(i (b) energy region for iron nitride 
Fe4N, e/{a.u.f 
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6 Summary 



We performed ab initio band structure calculations of iron nitride Fe4N. All our re- 
sults are in very good agreement with experimental data. Structure characteristics 
as well as magnetic properties are well reproduced. We also performed calculations 
of hypothetical Fe4N disordered structure. Comparative analysis of these two calcu- 
lations showed that the formation of Fe-N pairs in Fe4N compound is energetically 
unfavourable. 

We also calculated hyperfine interaction parameters for both Fe4N structures tak- 
ing into account the presence of magnetic field. For iron nitride structure we obtained 
perfect agreement of quadrupole splitting values with experimental data. For disor- 
dered Fe4N structure we observe considerable decreasing of quadrupole splitting for 
Fe-N-Fe-N chains, and very small quadrupole splitting for Fe-N pairs. 

Detailed analysis of EFG formation in iron nitride Fe4N allows us to conclude that 
EFG is mainly formed by space distribution of nitrogen 2s and 2p electrons, although 
it is greatly canceled by Fe3ci— electrons. Magnetization reduces this cancellation effect 
and increases the value of total EFG in iron nitride structure. 
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